Efficient Iterative Solvers for Stochastic Galerkin Discretizations of Log-Transformed Random Diffusion Problems
نویسندگان
چکیده
We consider the numerical solution of a steady-state diffusion problem where the diffusion coefficient is the exponent of a random field. The standard stochastic Galerkin formulation of this problem is computationally demanding because of the nonlinear structure of the uncertain component of it. We consider a reformulated version of this problem as a stochastic convection-diffusion problem with random convective velocity that depends linearly on a fixed number of independent truncated Gaussian random variables. The associated Galerkin matrix is nonsymmetric but sparse and allows for fast matrix-vector multiplications with optimal complexity. We construct and analyze two block-diagonal preconditioners for this Galerkin matrix for use with Krylov subspace methods such as the generalized minimal residual method. We test the efficiency of the proposed preconditioning approaches and compare the iterative solver performance for a model problem posed in both diffusion and convection-diffusion formulations.
منابع مشابه
Solving log-transformed random diffusion problems by stochastic Galerkin mixed finite element methods
Stochastic Galerkin finite element discretisations of PDEs with stochastically nonlinear coefficients lead to linear systems of equations with block dense matrices. In contrast, stochastic Galerkin finite element discretisations of PDEs with stochastically linear coefficients lead to linear systems of equations with block sparse matrices which are cheaper to manipulate and precondition in the f...
متن کاملEfficient Solvers for a Linear Stochastic Galerkin Mixed Formulation of Diffusion Problems with Random
We introduce a stochastic Galerkin mixed formulation of the steady-state diffusion equation and focus on the efficient iterative solution of the saddle-point systems obtained by combining standard finite element discretizations with two distinct types of stochastic basis functions. So-called mean-based preconditioners, based on fast solvers for scalar diffusion problems, are introduced for use ...
متن کاملEfficient Solvers for a Linear Stochastic Galerkin Mixed Formulation of Diffusion Problems with Random Data
Abstract. We introduce a stochastic Galerkin mixed formulation of the steady-state diffusion equation and focus on the efficient iterative solution of the saddle-point systems obtained by combining standard finite element discretisations with two distinct types of stochastic basis functions. So-called mean-based preconditioners, based on fast solvers for scalar diffusion problems, are introduce...
متن کاملIterative Solvers for the Stochastic Finite Element Method
This paper presents an overview and comparison of iterative solvers for linear stochastic partial differential equations (PDEs). A stochastic Galerkin finite element discretization is applied to transform the PDE into a coupled set of deterministic PDEs. Specialized solvers are required to solve the very high-dimensional systems that result after a finite element discretization of the resulting...
متن کاملA Kronecker Product Preconditioner for Stochastic Galerkin Finite Element Discretizations
The discretization of linear partial differential equations with random data by means of the stochastic Galerkin finite element method results in general in a large coupled linear of system of equations. Using the stochastic diffusion equation as a model problem, we introduce and study a symmetric positive definite Kronecker product preconditioner for the Galerkin matrix. We compare the popular...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- SIAM J. Scientific Computing
دوره 34 شماره
صفحات -
تاریخ انتشار 2012